On the low frequency electromagnetic waves in quasi-two-dimensional metals 
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Here, we theoretically analyze spectra of weakly attenuated electromagnetic waves which may 
appear in a Fermi-liquid of charge carriers in quasi-two-dimensional (Q2D) layered conductors when 
an external magnetic field is applied perpendicularly to the conducting layers. We study transverse 
modes propagating along the magnetic field. The frequencies of the modes are assumed to be lower 
than the cyclotron frequency of the charge carriers. It is shown that Fermi-liquid interaction of the 
charge carriers in Q2D conductors gives rise to a mode which cannot appear in a gase of charged 
quasiparticles, as well as it happens in conventional metals. Also, we show that the Fermi surface 
plofile may cause significant changes in the waves spectra and we analyze these changes. 

PACS numbers: 7l.l8.+y, 71.20-b, 72.55+s 
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In the last two decades various quasi-two-dimensional 
(Q2D) materials with metallic type conductivity such as 
organic metals, intercalated compounds and some other, 
have been synthesized. Galvanomagnetic phenomena 
and quantum oscillatory effects in these materials were 
intensely studied (see e.g. Refs. [l], 0, 0, HI)- It is 
a common knowledge that manifold weakly attenuating 
electromagnetic waves may appear in the electron liq- 
uid of conventional three dimensional metals in the pres- 
ence of a strong magnetic field @, Q. The conditions 
for the occurence of the electromagnetic waves may be 
much more favorable in Q2D conductors than in conven- 
tional metals. This happens due to a strong anisotropy 
of electrical conductivity in such materials. As recently 
shown [f| , the specifics of the Q2D electron spectrum may 
cause peculiar features in the dispersion of some collec- 
tive modes whose frequencies to are lower or of the order 
of the cyclotron frequency of the charge carriers fi. 

In the present work we mostly concentrate on the effect 
of the Q2D electron spectrum on the transverse Fermi- 
liquid (FL) cyclotron mode first predicted by Silin Q. 
The very existence of this mode results from the FL corre- 
lations of conduction electrons. In a metal with a nearly 
spherical Fermi surface (FS) this wave is circularly polar- 
ized, and it propagates along the applied magnetic field. 
We start our consideration by using the simple tight- 
binding approximation for the Q2D electron spectrum: 
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Here, z axis is assumed to be perpendicular to the con- 
ducting layers; pj_, p z are the electron quasimomentum 
projections on the layer plane and z -axis, respectively; 
the effective mass m corresponds to the electron motions 
in the layer plane; d is the interlayer distance. The pa- 
rameter £ in the Eq. (1) is the interlayer transfer inte- 
gral. We restrict our analysis with strong but nonquan- 
tizing magnetic fields Ti/t <C Tifl <T (r is the electron 
scattering time, and T is the temperature expressed in 
energy units). 



Within the phcnomcnological Fermi-liquid theory 
electron-electron interactions are represented by a self- 
consistent field affecting any single electron included in 
the electron liquid. Due to this field the electron energies 
E(p) get renormalized, and the renormalization correc- 
tions depend on the electron position r and time t : 



AE = Tr„ 



(27T?i) = 



F{p,a-p',a')5p{p',v',B',t). (2) 



Here, 6p(p, r, &,t) is the electron density matrix, p is 
the electron quasimomentum, and a are spin Pauli ma- 
trices. The trace is taken over spin numbers a. The 
Fermi-liquid kernel included in Eq. (2) is known to have 
a form: 

F(p, a; p', &') = p(p, p') + 4(&a')V(P> p') (3) 

For an axially symmetric FS the functions f{p,p') and 
ip(p, p') actually depend on the angle between pj_ 
and pi and on the longitudinal components Pz_,p' z - Ac- 
cordingly, we may present ip(p, p') as follows [9(: 

ViPziP'z, cos 0) = ipoo+PzP' z foi + (P±P±)( l Pio+PzP' z ( Pii)- 

(4) 

The function ip(p,p') may be presented in the similar 
way. In the Eq. (4) the functions ipoo, tpoi, <fio, fu are 
even in all their arguments, namely: p z ,p' z and cos 6. 

Also, we assume that both magnetic field B = (0, 0, B) 
and the wave vector of the considered wave q = (0, 0, q) 
are directed along the FS axis. Then the response of the 
Q2D electron liquid could be expressed in terms of the 
electron conductivity circular components a± (u>, q) = 
o~ xx (lj, q) ±«CT yx (o;,q). Substituting the expression (4) 
for the FL function </?(p,p') into the linearized transport 
equation for the nonequilibrium distribution function for 
the conduction electrons g(p, r,t) — Tr a (Sp(p,r,a,t)) 
one may derive the expressions for a± (cj, q) including 
the terms originatin g fro m the FL interactions. As shown 
in earlier works \v\. for an axially symmetrical FS 
and within the adopted geometry (both q and B being 
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in parallel with the FS symmetry axis) the conductivity 
components er±(u>,q) take the form: 



2me p F 



Here, 



4> d 



1 a{x)x n dx 
■ l U X± T v(x) ' 



(5) 

(6) 
(7) 

(8) 



<3n = / a(x)x n dx; 



X± = 1 ± fi/w + i/iOT] 
a(x) = A(x)/A Q ; 



x = Pz/po; 
v(x) = v z /v ; 



where vo,po are the maximum values of the longitudi- 
nal components of the charge carriers velocity and quasi- 
momentum, respectively: Aq = irpp; p = mvp — 
\f2mEp; Ep is the Fermi energy. The coefficients a>i : 2 
characterize FL interactions of the charge carriers: 



0-1,2 ■ 

A = 

h = 



A 



1 + A,a' 
2m r po 



2m 
(2nTi) 3 



p^tpwdpz 



/Po 
P±<pupldpz 
-po 



(9) 



When the charge carriers spectrum is described with 
the Eq. (1), the integrals Q n take on the form: 



1 [1 + ricos(-Kx)]x n dx 

-i U X± T sin(7ra;) 
l 

[1 + 1] cos(irx)]x n dx; 



(10) 
(11) 



where r) = 2(/Ep <C 1. Also, v = irrjVFpp /po <C vp. 

The dispersion equation for transverse modes propa- 
gating along the magnetic field could be writen as fol- 
lows: 
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c q 



Aftiwa-i-{<jj, q) = 0. 



When the relevant charge carriers are electrons, this 
equation has solutions corresponding to weakly attenuat- 
ing helicoidal and Fermi-liquid cyclotron modes for " — " 
polarization. Assuming the charge carriers to be holes, 



one must choose " + " polarization to get similar solu- 
tions. For certainty, we further analyze collective modes 
in the electron liquid. 

Using the expressions (10), (11) we may straightfor- 
wardly calculate the integrals $~ and Q n . For instance, 
we get 
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$0 = - ((«X-) 2 - 1) 



-1/2 



(13) 



Then, substituting the obtained results into the equa- 
tion (12) we find the solutions corresponding to the low 
frequency waves within the collisionless limit (r — > oo) : 



U>1 



Wg 2 -(g 2 -i)(/i-&7g 2 ) 2 
i - (A - W) 2 



Here. 
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a(x)v 2 (x)x 2 dx — —— 
l Qo 
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(14) 



(15) 



1 x 2 

a(x)v(x)xdx 



2 V 2ir 2 



(16) 
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q = qvo/fl; uj = w/Cl; and uj p is 

the electron plasma frequency. These solutions describe 
weakly attenuated collective modes in the Fermi-liquid 
of a Q2D conductor provided that the FS warping is not 
too small (r] > (gZ) _1 , where I = vpr is the electrons 
mean free path in the conduction layers planes). Within 
the limit 77 — ^ which corresponds to a 2D metal with 
a purely cylindrical FS, the expressions (14), (15) are 
reduced to <Di = 1 + fi , 0J2 = 1 + A, respectively, and 
they do not describe any electromagnetic waves. This 
is reasonable, for electromagnetic waves cannot travel in 
the electron liquid of a layered conductor perpendicularly 
to the layers when the interlayer transfer integral £ (see 
Eq. (1)) becomes zero. 

The solution (14) agrees with the corresponding result 
obtained by Kirichenko et al (see Ref. Q). It describes a 
helicoidal wave whose spectrum at small q (q <^ 1) has 
a form: 
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The spectrum of the wave is changed due to the FL inter- 
actions whose effect leads to the occurence of the term A 
in Eqs. (14), (17). In conventional 3D metals vq ~ vp , 
(12) therefore mvo/po ~ 1. So, assuming u 2 /c 2 q 2 ^> 1, we 
may arrive at the well known expression for the frequency 

13. The heli- 



of the helicoidal wave: to = Qc q /lo 2 
coidal wave may appear in an electron gas, and FL inter- 
actions manifest themself only through slight changes in 
the wave dispersion. 
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FIG. 1: Dispersions of low frequency transverse electromag- 
netic waves traveling along the magnetic field directed per- 
pendicularly to the conducting layers in a Q2D conductor. 
The diagonal straight line corresponds to the Doppler-shifted 
cyclotron resonance. Left panel: The helicoidal (lower curve) 
and the Fermi-liquid (upper curve) modes. The curves are 
plotted using Eqs. (14), (15) for r) = 0.075, /i = 0.2, f 2 = 
—0.2, B = IT. Right panel: The dispersion curves for heli- 
coidal waves at various values of the parameter r\. The curves 
are plotted for r/ = 0.025, 0.05, 0.075 (from the left to the 
right), at Qi = 0.2, B = IT. 
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FIG. 2: The effect of the FS shape on the dispersion of 
the Fermi-liquid cyclotron waves. The curves are plotted at 
ai = 0.2, a 2 = -0.2, B = IT. Left panel: Fermi-liquid 
cyclotron waves at various values of the parameter r\ for the 
cosine- like rippled FS corresponding to the energy-momentum 
relation (1). The curves are plotted for r\ = 0.1, 0.075, 0.05 
and 0.025 (from the left to the right). The dash-dot line 
indicates values of the cutoff wave vector q m for all disper- 
sion curves. Right panel: The helicoidal (lower curve) and the 
Fermi-liquid (upper curve) waves in a Q2D metal with a com- 
plex FS profile. The curves are plotted using Eq. (20) with 
s = 2 to describe v(x), at /i = 0.2, f 2 = -0.2, -q = 0.075. 



On the contrary, another solution of the dispersion 
equation describes a mode which exists due to the FL 
interactions, and cannot occur in the electron gas. The 
edge frequency of this mode is close to the cyclotron fre- 
quency: uj(0) = 0(1 + /2). As already mentioned, such 
Fermi-liquid cyclotron modes were predicted and their 
spectra analyzed by Silin for an isotropic electron liquid. 
Within the isotropic model vq = vp, pq — pp, /U = 8/35, 
so the expression (15) matchs the corresponding results 
of Ref. [8|| . For certainty, we further assume that the 
FL parameter f 2 takes on a negative value. Then the 
frequency of the Fermi-liquid cyclotron wave is smaller 
that the cyclotron frequency f2. 

Both solutions of the dispersion equation (12) are 
shown in the Fig. 1 (left panel). The diagonal straight 
line in the figure corresponds to the Doppler-shifted cy- 
clotron resonance uj — 1 — q. This line separates the 
region where the collisionless Landau damping mecha- 
nism works (the domain above the line) from the region 
where it does not work (the domain below the line), and 
where weakly attenuating electromagnetic waves may oc- 
cur. The curves in the figure are plotted according to Eqs. 
(14), (15) assuming /1 = 0.2, f 2 = -0.2, r\ = 0.075 for 
B = IT. As shown in the figure, the dispersion curve 
for the Fermi-liquid mode is located within the region 
free from the collisionless damping at small q. When q 
reachs the cutoff value q m ~ \a 2 \, the curve meets the de- 
marcation line uj — 1—q, and it terminates here. Beyond 
this line the wave is damped due to Landau collisionless 
attenuation. As for helicoidal mode, its frequency is sig- 
nificantly smaller than f2, and the wave vector may ac- 



cept any value within the range < q < 1. At q ~ 1 
the dispersion curve of this mode is located very close to 
the line Q = 1 — q which corresponds to the Doppler- 
shifted cyclotron resonance. Electron collisions give rise 
to the attenuation of the wave. However, the attenuation 
remains weak when (fir) -1 < q < 1 — (fir) -1 . 

The effect of the FS shape on the dispersion curves of 
helicoidal modes is shown in the right panel of the Fig. 
1. We see that the maximum frequency becomes greater, 
and it is reached at smaller values of q when the FS 
warping decreases. In the limit r/ — ^ the curve is to 
be reduced to the point q — 0; Q = 1 + f\. Also, the 
FS warping has a noticeable impact on the dispersion 
of the Fermi-liquid wave, as demonstrated in the Fig. 2 
(left panel). Here, the dispersion curves are plotted at 
f 2 = —0.2, and r\ takes on values from 0.025 to 0.1. 
The cutoff values of the wave vector magnitude q m corre- 
spond to the points of intersection between the dispersion 
curves and the dash-dot line in the figure. We see that 
for weaker FS warping the value of q m increases, so the 
wave may appear at greater wave vectors (shorter wave- 
lengths) provided that the magnetic field is held constant. 
To bette estimate the effect of a Q2D electron spectrum 
on the Fermi-liquid cyclotron mode, we may compare the 
results for the cutoff value q m presented in the Fig. 2 
with those obtained for quasi-isotropic 3D metals. It is 
shown that for such metals q m vp/il ~ I/2I (see |13|). 
Therefore, at f% = —0.2, q m vp/tt ~ 0.2 in a quasi- 
isotropic metal, whereas in Q2D metals with a moder- 
ately warped FS q m vp/Q may reach values significantly 
greater than unity at the same value of f%. 
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The energy- momentum relation (1) is a simplest ap- 
proximation used to describe electron spectrum in Q2D 
conductors. Within this model the maximum value of 
the longitudinal velocity is reached at p z — ±po/2 (x — 
±1/2) which correspond to the inflection lines on the FS. 
Near this lines v(x) may be approximates as: 

?(x) = ilTy( lT l/2) 2 . (18) 

In realistic Q2D metals the profiles of nearly cylindrical 
FSs may differ from a simple cosine warping provided 
by the Eq. (1). In general, neglecting the anisotropy in 
the layers planes, we may present the energy-momentum 
relation in the form [Bj]: 

n—l v / 

where e„ are dimensionless coefficients. By using this ex- 
pression we get opportunity to describe various FS pro- 
files, varying the values of e n . In particular, assuming 
t\ = 1 and all remaining e n = 0, we reduce the general 
expression (19) to the form (1). 

Assuming that the shape of the rippled FS cylinder 
may differ from that described by Eq. (1) we may con- 
jecture that in some Q2D conductors the longitudinal 
velocity v(x) near the inflection lines x — ±x* can be 
presented as follows: 

v(x)=±1tp 2 (xTx*) 2s (20) 

where s > 1, and p 2 is the dimensionless coefficient of 
the order of unity. The parameter s characterizes the 
shape of the FS segments where the longitudinal velocity 
v(x) gets maximum/minimum values. When s > 1 the 
relative number of the electrons moving with the max- 
imum longitudinal speed becomes significantly greater 
than in the case of a cosine-like FS profile (see eq. (18)). 
These electrons take a major part in the formation of the 
transverse collective modes near the Dopplcr-shifted cy- 
clotron resonance. The increase in their number brings 
changes in the spectra of such waves. The changes are 
especially significant for the Fermi-liquid mode. 

The effects of the FS local geometry on the spectra 
of the transverse Fermi-liquid waves traveling along the 
magnetic field in conventional 3D metals were studied in 
the earlier work Q . It was shown that when v(x) is de- 
scribed by the expression of the form similar to the Eq. 
(20), the dispersion curve for the Fermi- liquid mode does 
not terminate at q — q m - Instead, the curve has a con- 
tinuation in the region of shorter wavelengths where its 
frequency becomes low (u> <C 1). These conclusions re- 
main valid for Q2D conductors. When s takes on values 
greater than unity, the dispersion curve of the Fermi- 
liquid cyclotron wave runs below the demarcation line 
uj = 1 — q up to q = 1, as shown in the figure 2 (right 



panel). The low grequency part of the wave spectrum is 
located very close to this demarcation line. The latter 
corresponds to the Doppler-shifted cyclotron resonance, 
as it was mentioned above. 

Finally, in the present work we studied the propaga- 
tion of low frequency (u) < f2) electromagnetic waves in 
a Fermi-liquid of charge carriers in Q2D conductors. We 
analyzed transverse modes propagating along the mag- 
netic field applied along the normal to the conducting 
layers. We found that the mode arising due to FL in- 
teractions of charge carriers may occur in Q2D metals 
alongside the helicoidal mode analyzed in the Ref. [5J . It 
appears that this Fermi-liquid cyclotron mode (and the 
helicoidal mode, as well) may be easier observed in Q2D 
conductors than in usual 3D metals. This happens due 
to the smallness of the charge carriers velocity compo- 
nent along the normal to the conducting layers in Q2D 
metals. The smallness of v z causes the possibility for 
the collective modes to exist within a broader range of 
the wave vector and wavelength magnitudes than in usual 
3D metals. Also, we briefly analyzed the effects of the FS 
local geometry in Q2D conductors on the dispersion of 
the Fermi-liquid cyclotron wave, and we showed that for 
some FS profiles the mode may exist at frequencies which 
are significantly lower than the cyclotron frequency. 
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